Order and supersymmetry at high filling zero-energy states on the triangular lattice 
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We perform exact diagonalization studies in d — 2 dimensions for the Fendley and Schoutens 
model of hard-core and nearest-neighbor excluding fermions that displays an exact non-relativistic 
supersymmetry. Using clusters of all possible shapes up to 46 sites, we systematically study the be- 
havior of the ground state phase diagram as a function of filling. We focus on the highly degenerate 
zero-energy states found at fillings between 1/7 and ~ 1/5. At the lower end of that interval, at 
filling 1/7, we explicitly show that the ground states are gapped crystals. Consistent with previous 
suggestions, we find that the extensive entropy of zero states peaks at a filling of ~ 0.178. At the 
higher end of the interval, we find zero energy ground states at fillings above 1/5, in contradic- 
tion to previous numerical studies and analytical suggestions; these display non-trivial amplitude 
degeneracies. 



I. INTRODUCTION 

Supersymmetry is an extended symmetry which part- 
ners together fermions and bosons and is typically stud- 
ied in the context of extensions of the high-energy Stan- 
dard model of particle physics^. Recently, Fendley and 
Schoutens^l proposed a many-body lattice model which 
exhibits what is known as supersymmetric quantum me- 
chanic sSS. It consists of hard-core lattice fermions with 
nearest neighbor exclusion and a precisely fine-tuned in- 
teraction coupling. The eigenstates of the model consist 
of positive doubly degenerate energy states where the 
supersymmetry is spontaneously "broken" and unpaired 
zero-energy ground states. The nature of these states 
was studied extensively on the square lattice^, in spe- 
cially decorated two-dimensional latticed^, and also in 
one-dimensional models, where it provided insights into 
hidden many-body symmetries. 

On the square lattice the number of zero-energy states 
was exponential only in the system's linear dimension 6 , 
leading to a subextensive entropy of possible ground- 
states. In contrast, on the triangular lattice the num- 
ber of zero energy states appears exponential in the area, 
as recent exact diagonalization (ED) studies have indi- 
cated^. Jonsson, by studying finite triangular clusters, 
recently conjectured^, that on the triangular lattice, 
zero-energy states only appear in the interval 1/7 < / < 
1/5, where / is the filling (number of fermions/site)P3 
Huijse et aZP^ recently investigated, using ED, the na- 
ture and number of zero-energy states on the triangular 
lattice. Their results appear consistent with Jonsson's 
conjecture and the suggestion of exponential zero energy 
states' degeneracy, but they focused on quasi-one dimen- 
sional ladders, of width up to four. 

In this paper, we use exact diagonalization applied to 
all symmetrically inequivalent periodic clusters of vari- 
ous shapes and sizes, to determine the phase diagram as 
a function of filling. By combining data from many dif- 



ferent clusters, we are able to obtain curves of the energy 
as a function of filling across the whole phase diagram, 
and the entropy of zero states as a function of filling, 
across the interval in which they are found. Our study 
emphasizes too-dimensional clusters, wider than those of 
RefP. 

The outline of our paper and main findings are as fol- 
lows. First, in Section [TTJ we introduce the supersym- 
metric lattice model and the parameters of our exact di- 



agonalization calculations. In Section III we present our 
results on the different phases that appear as the density 
of fermions increases: a Fermi liquid phase at low filling 
/ < 1/7, zero-energy states with extensive entropy at in- 
termediate fillings, and high-density states (with tenden- 
cies to spatial orders) from /« 1/5 up to the maximum 
filling / = 1/3. In Section IV we focus on the states 



found near either endpoint of the zero-energy ground 
state interval. The minimum filling with a zero energy 
state is exactly / = 1/7 as conjecturecP} however, from 
that filling up to / w 0.156, the only zero-energy states 
are crystal-like and gapped states, and have no extensive 
entropy. At fillings slightly higher than 1/5, we find that 
certain zero-energy states exist, contradicting the con- 
jecture^-H These states show a tendency to anisotropic 
forms of spatial order, and their wavefunctions exhibit 
surprising regularities, in that many inequivalent fermion 
configurations have the same, maximal amplitude in the 
wavefunction (Sec. IVB3I. Finally, Section [vj contains 
a summary of our results and speculations on how these 
results may fit into a complete picture of the phase dia- 
gram. 



II. MODEL AND METHOD 



The Hamiltonian of the model is based on the defini- 
tion of the operators Q = J2 c \P(i) an d = J2 c iP(i)' 



where P, 



(i) 



n 



j next to i v 



j(l — CjCj) form projectors that 



exclude nearest neighbor occupancy. Following the ba- 
sic recipe on constructing the supersymmetric quantum 
mechanics^, the Hamiltonian is just 

H = {Q\Q} = tY J E p <*> c I^->+^E p <*>> 

i j next to i % 

(i) 

with t = V . The first term in ([I]) is a kinetic term of 
nearest-neighbor hopping (subject to the hard-core con- 
straint); the second term is a potential energy, equal to 
the number of fermions plus the number of vacant sites 
not forbidden by a neighboring fermion, which effectively 
includes pairwise farther-neighbor interactions as well as 
multi-fermion interactions. The Hamiltonian ([!]) is su- 
persymmetric only when finely tuned to the special point 
t = V. In view of the massive ground-state degeneracy , 
any change in the Hamiltonian would be a singular per- 
turbation leading to a new phase, if the energy spectrum 
is gapless; hence, the t = V point in the extended phase 
diagram is speculated to be a multicritical point at which 
various phases meelP. 

The model contains an exact Af — 2 supersymmetry 
in that all eigenstates with nonzero energy E belong to 
degenerate pairs, \a) with F fermions and \b) with F + 1 
fermions, such that Q\a) = y/~E\b) and Q^\b) — yE\a), 
whereas Q\b) = Q'\a) = 0. The state with an even num- 
ber is considered "bosonic" and the one with an odd num- 
ber is considered "fermionic" . An interesting corollary is 
that if there is a branch of fermionic elementary excita- 
tions, there must be a corresponding bosonic branch and 
vice versa. 

Most importantly, there may be ground states \G) of 
zero energy, such that Q\G) = Q'\G) = 0. A thermo- 
dynamically extensive zero-energy ground state entropy 
was reported and analyzed in several lattices, including 
integrable chains or ladders^^ as well as two-dimensional 
lattices, of which the triangle lattice is the simplest. In 
each model, the zero-energy states are always limited to 
fillings in the interval / G (/_, /+), where / is the num- 
ber of fermions per site. These prior studies mostly indi- 
cated that the zero-energy states have a strong tendency 
to form crystalline phases. 

We diagonalize the model of Eq. ([IJ on a triangular 
lattice using clusters of different size and shapes. Each 
cluster is characterized by the shortest and the next to 
shortest edge, v = (vi,v 2 ) and u = (ui,u 2 ) respectively, 
where the integer coordinates ui <2 and v\ t2 are expressed 
relative to the triangular lattice basis vectors e\ — x and 
e 2 = — |x+ the number of sites is N = u\v 2 — u 2 v\. 
We used all the symmetrically inequivalent clusters such 
that v belongs to the irreducible wedge of the triangular 
lattice and |v| 2 < |u| 2 for up to TV = 46 sites (larger for 
some fillings) . We perform the diagonalization separately 
for each sector defined by fermion count F and center of 
mass momentum. For the diagonalization we use the 
ARPACK packageP, which is an implementation of the 
Implicitly Restarted Lanczos Method (IRAM). Note that 
in accordance to Ref. [TUl well within the zero state in- 
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FIG. 1: Ground state energy per site vs. filling, minimized 
over all possible clusters with sizes 35 < N < A r max (/), where 
Nm^(f) = 56 for / < 1/7, 46 for 1/7 </ < 1/5 and 60 for 
/ > 1/5. To filter out ladders we also impose the constraint 
v > 4. For each filling the energy is minimized across all 
clusters and momenta. The different data points are colored 
according to the total system size. The solid line is a linear 
low envelop of the high energy regime which has the form 
E = | + A (/ - §), with A = 2.92. In the inset we show the 
energy versus filling for all the clusters zooming in the interval 
of zero energy states. 



terval we find that every momentum sector has a similar 
number of zero energy states, about 100 per sector in our 
largest system (N = 46); it is this huge degeneracy that 
limited the cluster size we could handle in this range of 
fillings. For other fillings, namely f < 1/7 and / > 1/5 
we study up to N — 56 and N = 60 respectively. In any 
case, for comparison, the tractable clusters are typically 
much larger than in Hubbard model exact diagonaliza- 
tion studies, since there is only one spin species and a 
hardcore constraint is enforced, each of which greatly re- 
duces the Hilbert space. 



III. RESULTS: PHASE DIAGRAM 

The model of Eq. |l]) displays three clear regimes as a 
function of particle density as evidenced from the equa- 
tion of state in Fig. [T] and in the configurations illustrated 
in Fig. [2] The qualitative behavior at low and high den- 
sity follows the studies by Henley and ZhangS^HSl w ho 
previously used exact diagonalization to study a model 
which was similar in that it has the same Hilbert space - 
spinless fermions with nearest-neighbor exclusion - but 
having V — in Eq. 0, so it lacked supersymmetry. At 
low densities a Fermi liquid phase appears, followed by 
zero energy ground states for fillings 1/7 </ < 0.2174. 
At the maximum allowed density of / — 1/3, the only al- 
lowed configuration is a crystalline phase with no quan- 
tum fluctuations, which, dominates the phase diagram 
beyond the interval of zero-energy states. In the rest of 
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FIG. 2: Idealized fermion configurations representing 
phases at various fillings. In this and subsequent figures, 
the triangular lattice sites are the centers of the hexagons, and 
filled hexagons are spinless fermions. For each filling, a config- 
uration of minimal potential energy is shown. More generally, 
depending on the unit cell constraints of each particular clus- 
ter, the configurations of highest weight and having minimum 
potential energy tend to be mixtures of triangles from the 
crystals shown here, (a) At low density the system is a dilute 
gas of weakly interacting fermions. (b) A triangular super- 
lattice crystal at / = 1/7 minimizes the potential energy for 
that filling, (c) A covering configuration at / = 1/6 (d) Two 
possible covering configurations for / = 1/4 (e) Anisotropic 
triangular lattice at / = 1/5 (f) At / = 1/3, there are only 
three allowed configurations in which one of 3 sublattices is 
fully occupied and no hopping is possible. 



FIG. 3: Zero energy states entropy. The zero state 
entropy per site, averaged over all the clusters with the same 
filling / for different subsets of clusters. We imposed the 
constraint jvj > 4 to filter out the ladders which have a zero 
energy state at lower fillings. The solid line corresponds to 
all cluster with TV > 34 the solid green to clusters with 35 < 
TV < 41 and the dot dashed line to 35 < TV < 46. There seems 
to be a maximum in the vicinity of / = 0.178. 



approach would be the i-matrix formalism of Ref. 1171 



Zero-energy states (1/7 < / < /+) 



this section we discuss those three different regimes of 
the phase diagram. 



A. Fermi liquid phase (/ < 1/7) 

The low density Fermi liquid is s imilar to the 
one studied by Henley and ZhangpEI] in their non- 
supersymmetric model. The ground state for each fill- 
ing, at low density, has the quantum numbers and is ap- 
proximated by the corresponding ground state for non- 
interacting fermions on the same cluster. A consequence 
of the positive hopping t > in our Hamiltonian is that 
the single fermion dispersion has two symmetry-related 
valleys, with minima at the Brillouin zone corners. For 
a given cluster, the sequence of energies as fermions are 
added one by one shows "shell filling effects" as expected 
in a Fermi liquid^ that is, in the dilute limit two suc- 
cessive fermions have a similar addition cost (placing one 
in each valley), then the next two are higher. Thus the 
addition energies show an even/odd alternation on top 
of an overall increasing trend, due to the repulsive (hard- 
core) interaction. If it were necessary to incorporate in- 
teraction effects in a systematic way, one relatively simple 



The degeneracy of the zero-energy states is best char- 
acterized by their zero state entropy^, which is exten- 
sive in the intermediate fillings. This fact suggests that 
at the infinite system size limit, the wave function incor- 
porates a mix of local patterns (especially fragments of 
the crystals in Figure [2j However, at finite size particu- 
lar periodic boundary conditions may be commensurate 
with only a subset of these; hence, different clusters may 
show contrasting behaviors, which would all be present 
in different patches of a typical infinite-size configuration. 

Let A (F;C) be the number of zero energy states for 
each given cluster C; this is one of the outputs of the 
ED. We conjecture that A (F;C) ~ exp(TV5 (/)), where 
TV(C) is the number of sites in the cluster and / = F/N, 
with an entropy function Sq(/) that is well-defined in the 
thermodynamic limit. We construct a numerical approx- 
imation to the entropy function by first computing an 
entropy per site for each cluster and then averaging over 
all clusters (of any size) with the identical filling per site, 
that is 

S (f) = (j^log[No(fN(Cy,C)]^ . (2) 

This averaging is meant to reduce the commensurability 
effects of different cluster shapes. 
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To better approximate the entropy of a two dimen- 
sional system we excluded ladder-shaped clusters by 
putting a lower limit for the width of v. We found that 
ladder-like clusters tend to have larger entropies per site 
than the more two-dimensional ones, and also this en- 
tropy tends to peak at a smaller filling /. Including such 
clusters would have incorrectly biased our estimates and 
can spuriously suggest that So(f) has another peak below 
/ = l/6. 

The result is shown in Fig. [3| where we can can observe 
four interesting numerical trends 

(1) There are no zero-energy states below 1/7, verifying 
recent expectations^ and in agreement with recent 
numerical exact diagonalization studies.^ 

(2) In the interval 1/7 < / < 0.156 almost all clusters 
either have no zero energy states or have unique 
zero energy states in isolated momentum sectors. 
(In some cases there are doubly degenerate zero 
energy states due to rotational symmetries.) That 
accounts for the small fluctuations of the entropy 
in this range, visible in Fig. [3] In addition, as can 
be seen in the inset of Fig. [T] the zero energy states 
for / < 0.156 are gapped. 

(3) the bulk of the zero energy states lie in the in- 
terval 0.156 < / < 1/5. In this interval the ma- 
jority of the clusters have degenerate zero energy 
states in every momentum sector. Let (/o, S'max) 
be the maximum point of So(/) function; our re- 
sults indicate that /o 0.178 and S ma x ~ 0.152. 
Van EerteiF^, using a transfer matrix, found a nu- 
merical bound on the Witten index, which cor- 
responds^ to S'max > 0.13, and also indicates 
fo ~ 0.18, with which our results are consistent. 

(4) Finally, at the upper end of the interval of zero 
states in the vicinity of / ~ 1/5, the slope of the 
entropy function Sq(/) appears to become much 
steeper. Thus we cannot rule out the possibility 
that the zero-energy state interval may end below 
or at / = 1/5 in the thermodynamic limit. 

We have also made a study of the energy gaps in this 
filling interval in clusters that have zero states, i.e. the 
smallest non-zero eigen-energy (from any momentum sec- 
tor). The main prior result about gaps in this model on 
the triangular lattice is that ladders with / = 1/6 are 
gaplesg^l 

As just noted, zero-energy states with / < 0.156 ap- 
pear to be gapped. In the filling interval (0.156, 1/5) we 
attempted to study the gaps by averaging over different 
clusters, in the same spirit as Eq. However, although 
we see clear trends, the behavior of the gaps depends on 
the clusters in an apparently irregular way, so the re- 
sults arc trustworthy only within an order of magnitude. 
Since the gaps are strongly (but not always monotoni- 
cally) decreasing with N, we only average within groups 
of clusters having the same (F, N). 



Within our two-dimensional clusters having typically 
~ 40 sites, there seemed to be particularly low gaps 
around two fillings / > 0.170 and / s» 0.190, with a 
maximum around / ~ 0.179. It appears, in fact, that 
these two minima represent a supersymmetric pair with 
F and F+l fcrmions, e.g. fillings 7/42 and 8/42, and 
we conjecture that in the thermodynamic limit there is 
one minimum occuring at / ~ 0.178, where the entropy 
function is highest. The gap values at the minimum de- 
cay at least as fast as 1/2 F (up to F — 8 where our 
data is complete), and are 10 _3 -10~ 1 in our largest sys- 
tems. At other fillings around the middle of the inter- 
val (0.156, 1/5) we still see a decaying trend with F, but 
slower. For / < 0.163 and / > 0.195 gaps larger than 
0.1 are seen even in our largest systems, so we cannot 
definitely say whether these are gapped. 



C. Phase diagram at filling / < 1/3 

It is gencrically expected that strongly interacting lat- 
tice fcrmion models at high densities phase separate into 
a high-density insulator, such as the inert / = 1/3 crystal 
in Figure [2](f ) , and a low-density liquid phase in which 
kinetic energy is dominant. On the triangular lattice, 
we would guess that the coexisting liquid density falls 
around / = 0.21 where the Hilbert space is largest. (We 
used the Pauling approximation for the entropy of al- 
lowed configurations, following the Appendix of Ref.HU) 
That is roughly the filling at which the dense liquid is be- 
coming congested, hence less favorable, due to the hard- 
core constraint. 

Just which dense st ate d oes this dense liquid coex- 
ist with? Prior researc h 1 14 ! 15 ! on the non-supersymmetric 
model on the square lattice showed that the dense phase 
is not the maximally occupied crystal. Instead, it is a 
crystal containing a dilute array of quantum-fluctuating 
strings, which we call "stripe-walls" since each has a 
deficit of fcrmion density and is a domain wall of the crys- 
tal ordeJ^l There is no previous literature on such stripe- 
walls in the case of the triangular lattice for fermions, 
only for hardcore bosons ( in th e context of 4 He ad- 
sorbed on carbon nanotube d 2 ° l 21 l Nevertheless, it was 
shown in Ref. [T3] that in an isolated stripe-wall, hard- 
core constraints do not allow particle exchanges, hence 
the fermion and boson cases should be in fact equivalent. 

The stripe- wall runs perpendicular to one of the bond 
directions; fermions on the edge of a domain are free 
to hop in that direction, owing to the relative shift of 
the other domain. The deficit of electrons is 1/3 per 
step along the stripe-wall. Let the energy per step be 
-Esw : then the chemical potential associated with (non- 
interacting) stripe- walls is nsw — 3i?sw- Thus, a phase 
of dilute stripe walls is represented by an energy function 
E(f) = 1/3 - msw(1/3 - /)• #(l/3)=l/3 in the maxi- 
mally filled crystal, since there is no kinetic energy and 
the potential energy in ([I]) is F = N/3. 

Note that as / continues to decrease and the array 
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becomes less dilute, collisions between adjacent stripe- 
walls (due to their quantum fluctuations) start to become 
significant. This will typically reduce the kinetic energy, 
causing the E(f) curve to bend upwards. Thus, if /zsw 
exceeds the slope of a coexistence curve connecting to the 
/ = 1/3 crystal, the stripe- wall array is stable at / < 1/3. 
Its stability ends at the point where a tangent to the E(f) 
curve can be drawn to the coexisting liquid phase. (Such 
a point should exist since E{f) is upwards curving.) By 
contrast, if /^sw is less than the slope of the coexistence 
curve, the / = 1/3 crystalline phase coexists directly 
with a hole-rich phase having / s» 0.22. In the latter 
case, in the grand canonical-like ensemble, a first-order 
transition would be seen between the hole-rich phase and 
the / = 1/3 crystal. To check which scenario holds in 
the supersymmetric model, we must calculate the energy 
Esw- _ 

Past studies^22l showed that a single stripe- wall can be 
mapped exactly to a one-dimensional chain with nonin- 
teracting spinless fermions at half filling, with their hop- 
ping amplitude equal to t = 1. This mapping remains 
valid in the supersymmetric case, since all the accessible 
configurations for the stripe-wall have equal (and max- 
imal) potential energy F. Thus, the kinetic energy per 
step is — 2/tt, or —6/n per removed fermion; the poten- 
tial energy is +1 per removed fermion. Hence, we obtain 
/isw = (1 + 6/tt) = 2.91. 

For comparison, our numerics from exact diagonaliza- 
tion of the supersymmetric model showed that in the 
interval below / = 1/3, the energy is linear as a func- 
tion of filling with a form E(f) = § + A (/ — 4), with 
A = 2.92, represented by the solid line (cf. Fig. flj). Since 
we find that A = /isw, this leaves undecided whether a 
stripe- wall array is stable or the / = 1/3 crystal coexists 
directly with / m 0.22. 

The straight line fitting E(f) would pass through / ss 
0.22, close to the maximum filling of special zero-energy 
states we found above / = 0.2. However, these states do 
not appear in Figure [T] which is limited to robustly two- 
dimensional clusters; the numerical E(f) curve curves 
upwards in the vicinity of / = 0.22 and appears to hit 
the E = axis at / = 0.20-0.205. That means the zero- 
energy states have a smooth transition to E > states, 
which are stable in a short interval up to / = 0.22, and 
then possibly coexisting with either a stripe-array or the 
/ = 1/3 crystal. 



IV. ZERO-ENERGY STATES AT MAXIMUM 
AND MINIMUM FILLINGS 



The entropy of zero states vanishes, according to Fig- 
ure [3j at filling 1/7 and slightly over filling 1/5. Around 
those fillings, if zero states exist at all, they tend to be 
unique in their momentum sector and have other special 
properties, including a tendency to crystal-like spatial 
orderings. 
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56 
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(11,5) 


(2,6) 


56 
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(12,4) 


(1,5) 


56 
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(20,4) 


(1,3) 


56 


8 


(8,5) 


(0,7) 



TABLE I: The clusters that have a state at filling / = 1/7 
for 28 < N < 56. All the clusters are commensurate to the 
y/7 x \fl crystalline structure shown in Fig.J^b). In all cases 
there is a unique zero-energy state in 7 different momentum 
sectors with a leading configuration which resembles Fig.[2^b) . 




FIG. 4: A compilation of clusters with zero-energy 
states at filling / > \ for (a) v = (1,4), (b) v = (3,7) and 
(c) v = (2,4). For each case, we mark with different colors 
the vectors u of all the clusters with a high filling zero-energy 
state. The symmetrically equivalent u's use the same color 
(color online). The + and — signs correspond to center of 
mass momentum K = (0, 0) and = (0, n) respectively. 



A. Zero-energy states at filling / ~ 1/7 

The zero energy states near the lower end of the inter- 
val [/_,/+] have been discussed in Ref.UUlbut apparently 
not in detail. The authors, numerically, identified only 
one cluster, (7,0) x (0,7), having a zero energy state at 
/ = 1/7. In contrast, we found striking behaviors exactly 
at / = 1/7 for a variety of clusters. 

In Table [I] we present all the clusters with zero-energy 
states at / = 1/7. Interestingly, those are all of the 
distinct clusters which are commensurate with the \/7 x 
\fl crystalline order shown in Fig. [2|b), i.e. both the 
vectors u and v are lattice vectors of that pattern. The 
leading configurations of that wave function consist of the 
seven possible translations of the same perfectly ordered 
pattern. 

We find that those zero-energy states are gapped. 
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FIG. 5: Nature of the zero-energy states listed in Table |T] 
(a) The gap of the unique zero-energy state at / = 1/7 as a 
function of the number of particles F. The 3-leg ladders with 
v = (1,3) approach a different value for large F than the 
rest of the clusters. The geometry of the ladders is such that 
each unit cell of the v7 x V7 crystalline order wraps around 
it self, (b) The quantity q = 1 - (7|VWx| 2 ) 1/J? , the inferred 
probability per crystal site to find its fermion one site away, is 
independent of system size, as expected for long-ranged crys- 
talline order, (c) The weight of a sub-leading configuration 
(in which one fermion is off its site) relative to a leading one. 
This is also independent of F and approximately equal to 
hq/(l — q), as expected. 



Shown in Fig. 5(a) is the gap plotted against the clus- 
ter size, grouped as strips of different widths. We notice 
that the strip of shortest width, v = (1,3) [i.e. three-leg 
ladders], seems to have a gap approaching 0.46, whereas 
the strips with longer v have gaps larger than 1. We 
remark that the three-leg ladders are not representative 
of the thermodynamic limit, since the nearest neighbors 
transverse to the strip are images of the same fermion. 

It is possible to understand the crystalline order within 
a very simple picture, based on an eigenstate in which 
the symmetry is already broken. Our main hypothesis is 
that, to a good approximation, the probability \ip\ 2 , as a 
function of configuration, is a direct product of indepen- 
dent wavcfunctions of each fermion hopping in a separate 
hexagon of 7 sites, centered on one of the ideal crystal 
sites. Within each hexagon the weight is (1 — q) on the 
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FIG. 6: The leading configurations for two particular clus- 
ters in the interval 1/7 </ < 0.156. Both correspond to 
unique and gapped zero energy states at a particular momen- 
tum sector. Both of them have the form of stripes of the 
/ = 1/7 crystal separated by domain walls. 



center site and q/6 on the other sites. The leading con- 
figuration is when every fermion is on a center site and 
the leading weight is |?/>max| 2 = (1 — l) F Hi where the 
factor of 7 signifies that the actual eigenstate in any mo- 
mentum sector is a superposition of seven shifted copies 
of that eigenstate, w hich have no configurations in com- 
mon. Indeed in Fig. 5(b) we plot q = 1 — (7|V>max| 2 ) ^ 
and we find it independent of F with a limiting value of 
q « 0.148. 

A sub-leading term in the many-fermion wavefunc- 
tion (as we verified by inspecting the wavefunction) goes 
with a configuration in which each of (F — 1) fermions 
is on the center site of its hexagon, but one fermion 
is on a neighbor site in its hexagon. Hence the sub- 
leading amplitude should be IV'subl 2 = |f (1 — l) F ~ 1 an d 
IV'sublVlVwl 2 =q/6(l-q). According to Fig. [5(c) 



the 

limiting value of this ratio is 0.028 which corresponds to 
q rj 0.144. The good agreement of this to the q value 
mentioned above supports the validity of our simple pic- 
ture. 

The character of the / = 1/7 zero energy ground states 
survives in a small interval, 1/7 < / < 0.156 in which 
only some of the clusters have zero energy states, in only 
a few momentum sectors. Similarly to the / = 1/7 case, 
those states are found to be unique and protected by a 
gap which gets smaller as the / of these states gets larger. 
The inset of Fig. [T] shows that the energy vs filling does 
exhibit a gap in this interval. Looking at the structure of 
the wave function reveals that the leading configuration 
(and the symmetrically equivalent ones) constitute an 
f — 1/7 pattern with a domain wall as shown in Fig. [6] 
whereas the subleading configurations are those with one 
fermion hopped out of this pattern. 



B. Zero-energy states at filling / > 1/5 

The existence of zero-energy states for filling / > 1/5 
has been a subject of debate. Jonsson in a recent studjW 
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N 


F 


/ 


u 


V 


K 




VKl ea d 


AE 


28 


6 


0.2143 


(8,4) 


(1,4) 


(0,7r) 


9 


0.6667 


1.1038 


33 


7 


0.2121 


(9,3) 


(1,4) 


(0,0) 


12 


0.4893 


0.9327 


38 


8 


0.2105 


(10,2) 


(1,4) 


(0,7r) 


16 


0.3590 


0.8290 


43 


9 


0.2093 


(12,5) 


(1,4) 


(0,0) 


20 


0.2475 


0.7068 


48 


10 


0.2083 


(13,4) 


(1,4) 


(0,tt) 


25 


0.1705 


0.5978 


53 


11 


0.2075 


(14,3) 


(1,4) 


(0,0) 


30 


0.1127 


0.5075 


58 


12 


0.2069 


(16,6) 


(1,4) 


(0,tt) 


36 


0.0745 


0.4364 


29 


6 


0.2069 


(8,3) 


(1,4) 


(0,0) 


2 


0.3252 


0.2407 


34 


7 


0.2059 


(9,2) 


(1,4) 


(0,0) 


16 


0.4028 


0.6079 


44 


9 


0.2093 


(12,4) 


(1,4) 


(0,0) 


40 


0.2127 


0.3289 


54 


11 


0.2037 


(15,6) 


(1,4) 


(0,0) 


80 


0.1001 


0.3097 


44 


9 


0.2093 


(11,0) 


(2,4) 


(0,0) 


8 


0.1337 


0.7345 


44 


9 


0.2093 


(12,2) 


(2,4) 


(0,7r) 


8 


0.1337 


0.6338 


24 


5 


0.2083 


(6,6) 


(3,7) 


(0,0) 


4 


0.6379 


0.9461 


34 


7 


0.2059 


(7,5) 


(3,7) 


(0,0) 


14 


0.6198 


0.6386 


44 


9 


0.2093 


(8,4) 


(3,7) 


(0,0) 


32 


0.4144 


0.5823 


54 


11 


0.2037 


(9,3) 


(3,7) 


(0,0) 


66 


0.2531 


0.7621 



TABLE II: Clusters in which a high filling zero-energy state 
is found. iV is the number of sites, F the number of fermions, 
and / = F/N is the filling; the vectors u and v define the 
cluster unit cell. In each case, there is a unique zero-energy 
eigenfunction in momentum sector K. Aiead is the number 
of terms in that (normalized) wavefunction having identical 
maximum amplitude IV'icadl and not equivalent by transla- 
tions; Wi C ad = MeadlV'ieadl 2 is the combined weight of such 
terms; it can be seen that Wiead tends to vanish with increas- 
ing system size. Also, AE is the gap to the next state in the 
same sector (of F and K) 



in Figure|4j these clusters fall into sequences, each having 
fillings that approach / = 1/5 from above in the thermo- 
dynamic limit. 

The clusters with v = (1,4) in the beginning rows of 
Table IT] fall in a sequence with u = (0, 2 — e) + m (1, —1) 
modulo v for e = 0, 1 and m = 6, 7, . . . , , 12. These have 
F = m fermions in N = 5m + e — 2 sites so the filling is 
/ = /[5 + (e — 2)/m]. The momentum where the zero- 
energy states appear is K = (0, 0), except when e = and 
m is even we have K = (0,7r). (The e = 1 zero-energy 
states appear only for m = 6,7,9,11, so they may be 
a finite size effect.) Also, the last rows of the table are 
four clusters having v = (3, 7) and the other vector in a 
sequence u = 4 (1, 2) + m (1, —1) with m = 2, . . .. (the 
m = 1 cluster was omitted from the table because N < 
25.) These have F = 2m + 1 fermions in N — 10m — 1 
sites, for a filling which is obviously greater than 1/5. 

The gaps of the v = (1,4) ladders shown in table IT] 
can be fitted to a decaying exponential function of the 
system size, which extrapolates to a gapless state in the 
thermodynamic limit. Indeed, within this series th e gaps 
fit well to exp(— ConstA^ 2 ). We know (see Section IV A) 



that in a crystal the amplitudes should decay exponen- 
tially, so this is suggestive of a liquid phase, or perhaps a 
crystal with unbounded fluctuations, similar to a classical 
crystal in d = 2 at nonzero temperature. 

The clusters at width v = (3,7), the only ones in 
this table that are plausibly two-dimensional, have the 
largest gaps - of order unity - suggesting the presence of 
a gapped state in this filling. For the v = (2,4) ladders 
the data are insufficient to draw a conclusion. 



analyzed a collection of homology cycles on the triangular 
grid and showed that for a specific class of configurations 
there is an upper bound on the existence of zero-energy 
ground states which is / = 1/5. Moreover, Huijse et aZ™ 
through extensive exact diagonalization studies and ana- 
lytical arguments found agreement of the Jonsson conjec- 
ture with numerical results on the triangular grid. Our 
studies have found zero-energy states beyond the filling 
/ = 1/5. 

1. Occurrence and excitation gaps 

Figure [3] displays all the clusters where / > 1/5 zero- 
energy states are observed. Such states are observed in 
both elongated and roughly isotropic cluster shapes, sug- 
gesting that their existence is not an artifact of the aspect 
ratio (but it could be finite-size related). We find that 
for a few clusters there are unique zero-energy states at 
a single momentum sector which is either (0, 0) or (0, 7r) 
depending on the cluster. Table [TT| shows all the clusters 
and momentum sectors with / > 1/5 zero-energy states. 

The shapes of cluster that tend to support zero-energy 
states with / > 1/5 are shown in Table in] As illustrated 



2. Density correlations and coincident wavefunction 
amplitudes 

To further understand the ground state wavefunctions 
at fillings / > 1/5, we examined the density corre- 
lation functions in these wavefunctions [Fig. [7]ja,c,e)]. 
The density-density correlation function in each of these 
states suggests a non-trivial order, which in some cases it 
is stripe-like, whereas in others it corresponds to a more 
isotropic structure [cf. Fig. [7](a,c,e).] 

We also examined the leading amplitude configurations 
in each wavefunction [examples in Fig. [7][b,d,f)]. The 
combined weight of those leading configurations appears 
to vanish with the system size, but that happens even 
when the leading configurations strongly dominate the 
behavior, e.g. in the f — 1/7 crystal [see Figure [5^b)] 
where they decay as exp(— const N. For the / > 1/5 zero- 
energy wavefunctions, the decay appears to be faster, 
perhaps as s exp(— const A^ 2 . (This can be inferred from 
the last column of Table [TTJ ) This indicates that the 
/ > 1/5 states are less ordered than the / = 1/7 crystal. 




(e) (f) 

FIG. 7: The density density correlation (left panels) and 
the leading configuration (right panels) for the unique zero- 
energy states at / > 1/5. In each case, the unit cells are 
marked by light dotted lines, starting from the left lower cor- 
ner, note in (c,d) the second set of (1, 4) edges does not pass 
through the panel, while in (e,f) the second set of (2, 4) edges 
is the panel's left boundary, (a), (b) (8, 4) , (3, 7) at f = 9/44, 
(c), (d) (15, 6), (1,4) at / = 11/54, (e), (f) (12, 2), (2, 4) at 
/ = 9/44. (a), (c) and (d): map of the density density cor- 
relation (p(r)p(r + R)} as a function R. In each panel many 
copies of the unit cell (marked by thin lines) are shown; darker 
grays indicate higher correlations, except that the origin is 
made white for better contrast. At certain separations the 
correlation function is almost zero, but its only exact zeroes 
are the nearest neighbors to the origin (and points equivalent 
by translation), which are excluded by the reference fermion. 
(b), (d) and (f): one of the leading configurations all having 
identical weights as shown in table. [II] Each configuration has 
the covering property, i.e. every unoccupied site is neighbor 
to at least one fermion, and consequently the potential energy 
has a the minimum value. 



3. "Magtc" wavefunctions 

All of the / > 1/5 zero-energy wavefunctions have the 
unusual property of "coincident amplitudes", meaning 
there are multiple configurations having the identical am- 
plitudes (apart from sign). We have checked this only for 
those of maximum amplitude, i.e. the "leading" configu- 
rations. In all cases we know, these leading configurations 
all have the minimum possible potential energy V m i n . 

The coincident-amplitude property is somewhat rem- 
iniscent of the Rokhsar-Kivelson (RK) eigenstates^, in 
which all configurations have equa l weight. Generaliza- 
tions of the RK constructiorP^Hini might inspire future 
possibilities of analytic approaches to the supersymmet- 
ric model: e.g. one might construct the exact ground 
state wavefunction for the given cluster; or possibly, the 
quantum ground state might have the same weights as 
the Boltzmann ensemble of some classical model, which 
would assist in understanding the degree and nature of 
order in the thermodynamic limit. 

In the present work, we could only recognize the 
coincident-amplitude property when the zero state is 
unique/non-degenerate in its momentum sector; other- 
wise, the condition is ill-defined, since the amplitudes 
depend on which linear combination of the degenerate 
states is provided by the Lanczos solver. Thus, such 
wavefunctions are likely to be identified only near the 
endpoints of the interval of fillings with zero energy 
states, i.e. around / w 0.15 or / « 0.20; in practice 
we only identified them for / > 0.20. 

Furthermore, the leading configurations in Fig. (7^b) 
and (d) are dominated by fragments of the / = 1/5 crys- 
tal from Fig. [2^e), which is the densest simple crystal 
in which the fermions have available hoppings and the 
potential energy has the minimum value of F. In the ac- 
tual leading configurations, some triangles from the other 
crystals in Fig. [2] get mixed in, which allows more quan- 
tum fluctuations and lowers the kinetic energy. Notice 
these are very anisotropic structures, in that almost all 
allowed hoppings are in the same directions (aligned with 
the long direction of the / = 1/5 triangles). The config- 
uration in Fig. [r]^f ) appears to be dominated instead by 
rows of the / = 1/6 cell from Fig. [2jd) which explains 
why the (12,2), (2,4) cluster is more isotropic than the 
other clusters. 



V. CONCLUDING REMARKS 

The exact diagonalization study we performed on a va- 
riety of finite clusters illuminated certain aspects of the 
phase diagram of the supersymmetric lattice model of 
Eq. [I] The main findings are the presence of gapped 
crystalline zero energy states for / ~ 1/7, gapless quasi- 
ordered states for / > 1/5 in violation of Jonsson's con- 
jecture and a zero state entropy which peaks at / ~ 
0.178. Our results are not conclusive about the exact 
phase diagram in several places and we now speculate on 
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the possible scenarios, in light of our results. 

Our finding of an energy gap at / = 1/7, as shown in 
Section [IV A| implies that the Fermi liquid phase cannot 
connect continuously to the zero state at / = 1/7. Either 
there is a first-order transition, so that coexistence of the 
Fermi liquid with the / = 1/7 crystal is present in a small 
interval below / = 1/7, or else the crystal phase extends 
(with E > 0) to lower fillings. A specific mechanism 
for the latter alternative is that the crystal is doped with 
fermion vacancies (which thus form a Fermi liquid within 
the crystal), which destabilizes (i.e. melts) the crystal at 
a sufficient density of vacancies. Additional studies will 
be needed to decide between these scenarios. 

The next open question is the behavior immediately 
above / = 1/7. Our results suggest there is a family of 
states made by combining strips of the / = 1/7 and / = 
1/6 crystals: each such state is a zero-energy state similar 
to the / = 1/7 crystal in being gapped and in having 
a crystal-like symmetry breaking. This would imply a 
dense set of rational fillings above / = 1/7 that are, in 
some sense, isolated from each other rather than forming 
a continuum. 

We conjecture that the fillings in the interval 0.156 < 
/ < 0.2 are all gapless in the thermodynamic limit, based 
on the numerical observations discussed at the end of sec- 
tion |III B| If this is true in the thermodynamic limit, it 
would fit naturally with a liquid-like phase, in which the 
wavefunction has a vast number of important configura- 
tions, each of them mixing the units from all the crystals 
shown in Figure pi and the important configurations are 
accessible to each other by small steps corresponding to 
local rearrangements of these units. This scenario could 
be reconciled with the apparent increase of the gaps away 
from the center of that interval at / = 0.178, in the fol- 
lowing speculative picture. Say that the entropy (per 
site) of Zow-energy states has a similar dependence on / 
as the entropy of zero-energy states. Then, for / near the 
endpoints and with a finite size N, the expected number 
of low-energy states in a cluster might be less than unity, 
meaning that sometimes the cluster has none. 

For the clusters with zero energy states at / > 1/5 
we find non-degenerate zero-energy states in just a single 
momentum sector. The filling of those states seems to 
converge to 1/5 for increasing system size; furthermore, 
we checked that if we make a new cluster that simply dou- 
bles one of these clusters, and keep the same /, the larger 
cluster does not generically have a zero energy state. It 
is unclear how our results relate to the thermodynamic 
limit. These states showed more than one kind of order- 
ing tendency, depending on the cluster (as noted at end 
of Sec. IV B |, but the tendency is that most of the avail- 



able hoppings lie along a particular axis. This might be 
explained by a local order incorporating many fragments 
of the / = 1/5 crystal [Figure |2jje)], which only allows 
hops along one lattice axis. 

The corresponding wave functions at / > 0.2 are dom- 
inated by many configurations with equal ( "coincident" ) 



amplitudes that always belong to the lowest potential en- 
ergy sector, (Sec. IV B 3[ ). We certainly cannot rule out 
the possibility that similar wavefunctions exist for other 
zero states, since we are only able to detect this prop- 
erty when there is a unique zero state in a particular 
sector. The regularities in these wavefunctions suggest 
the possibility that some analytic structure, even an ex- 
act solution, may be found for the zero-energy ground 
states. 

Finally, at both ends of the range (1/5,1/3), we naively 
expect some form of coexistence of a hole-rich phase with 
the / = 1/3 inert crystal, or this crystal doped by an 
array of stripe- walls. The E{f) curve shows a straight 
line at very nearly the slope expected from a stripe- wall 
array. However, the picture from our exact diagonaliza- 
tions is not quite consistent with any scenario: the E(f) 
curve due to an array ought to curve upwards rather 
than be a straight line; in contrast, in the case of coexis- 
tence a straight tie-line is not expected in a finite system 
due to the additional cost of the domain wall. One in- 
triguing possibility is that the domain-wall cost is zero: 
that is, the phase at / > 0.5 and the stripe-wall array 
dissolve into each other. This is not implausible: the 
stripe- wall array is an / = 1/3 crystal with some units 
of the / = 1/5 crystal appearing along the stripe- walls, 
and it has hoppings entirely along one of the crystal axis 
directions. Thus, its properties are similar to some of 
the zero states we observed at / > 1/5, so perhaps the 
quasi-order we found in those states is simply a general- 
ization of the stripe- wall kind of state. It is worth noting 
that in Ref. [Till it was argued that the ground state at 
/ = 1/4 is gapped. It would be interesting to identify the 
deep origin of this gap, given that stripe- walls, as well as 
coexisting stripe-walls should be by all odds gapless. 

Overall, a possible interpretation of our consistent find- 
ing of a strong bias towards ordering in finite clusters, is 
that the supersymmetric model is a multicritical point 
in the parameter space of Hamiltonians, so that many 
different kinds of order may be stabilized by infinitesi- 
mal perturbations. The various sizes and shapes of the 
system introduce similar biases, so that different large 
clusters may be occupied by fragments of competing or- 
dered states which are all valid and degenerate in the 
thermodynamic limit. In future studies, it may be worth 
to add terms to the Hamiltonian ([I]) so as to explore the 
possible relation of supersymmetry to multicriticality^U. 
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